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ABSTRACT 

Experimental  observations  of  various  kinds  of  cavitation  instabilities  are  presented.  Typical  cavitation 
instabilities  occur  in  the  range  of  cavitation  number  where  the  inducer  head  is  not  affected  by  cavitation. 
This  type  of  cavitation  instability  is  caused  by  the  cavity  volume  increase  due  to  the  increase  of  incidence 
angle  and  is  treated  in  section  1.  Another  type  of  cavitation  instability  is  caused  by  the  positive  slope  of 
the  pressure  performance  curve  due  to  the  blockage  effect  of  cavitation.  This  type  of  cavitation  instability 
is  treated  in  section  2.  Stability  analyses  of  1-  and  2-dimensional  cavitating  flow  predict  various  modes  of 
cavitation  instabilities.  Usually  only  cavitation  surge  and  forward  rotating  cavitation  are  observed. 
Other  modes  of  cavitation  instabilities  such  as  backward  rotating  cavitation,  higher  order  modes  of 
cavitation  surge  and  rotating  cavitation  are  observed  only  under  limited  conditions.  Those  modes  of 
cavitation  instabilities  are  presented  in  section  3. 


1.0  TYPICAL  CAVITATION  INSTABILITIES  IN  A  3-BLADED  INDUCER 
1.1  Experimental  Apparatus 

The  pump  loop  used  is  shown  in  Fig.l  [1].  The  baseline  experiments  were  made  with  an  inlet  pipe  A  with 
an  inner  diameter  200mm  and  without  the  outlet  accumulator  B.  The  system  parameter  was  changed  by 
replacing  the  inlet  pipe  A  with  one  with  a  diameter  150mm  or  by  adding  an  accumulator  B  with  a  volume 

of  airl.5xl0-3  m3  under  the  test  condition.  The  base  pressure  (and  hence  the  cavitation  number)  was 
adjusted  by  using  a  vacuum  pump  connected  to  the  pressure  control  tank  with  about  0.65  m3  of  gas/vapor. 


Figure  1:  Test  Loop. 
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Figure  2  shows  the  test  section  and  the  performance  curve.  The  impeller  is  a  scale  model  of  LOX 
turbopump  inducer  for  HII  rocket,  with  three  blades  Zt  =  3 ,  outer  diameter  149.8  mm,  blade  angles  at  the 
tip  7.5°  at  the  inlet  and  9.0°  at  the  outlet;  the  housing  is  made  of  transparent  acrylic  resin,  with  inner 
diameter  150.8  mm  (constant  tip  clearance  of  0.5mm).  The  inlet  pressure  fluctuation  was  measured  at 
27.5mm  upstream  of  the  inducer  leading  edge,  at  two  circumferential  locations  with  various  separation 
angles  6.  The  number  of  cells  was  determined  from  a  plot  of  the  phase  difference  of  those  pressure 
signals.  The  rotational  frequency  of  the  impeller  was  maintained  N  =  4,000  ±  2 rpm  .  The  static  pressure 
performance  is  also  shown  in  Fig.2.  The  performance  curve  has  a  negative  slope  throughout  the  flow 
range  and  no  conventional  surge  nor  rotating  stall  is  expected  to  occur.  The  design  point  of  the  inducer  is 
(0  =  0.078,^=0.130).  Since  the  flow  meter  was  not  available  at  the  time  of  unsteady  pressure 
measurements,  the  operating  condition  is  shown  by  using  the  value  of  the  static  pressure  coefficient^ . 
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Figure  2:  Test  Section  and  Static  Pressure  Performance. 


1.2  Cell  Number  Identification 

Figure  3(a)  shows  an  example  of  cascade  plot  of  inlet  pressure  spectra  for  various  cavitation  numbers  o  , 
for  a  static  pressure  coefficient  y/s  =0.123  and  a  rotational  frequency  fN  =  4000/60  «  66.1Hz  .  Hence  the 
blade  passing  frequency  is  3  fN  =  200 Hz  .  We  observe  the  frequency  components  i  -  v  and  i'-v' .  If  we 
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represent  the  frequencies  of  i  and  V  components  by  f)  and  fv ,  we  observe  the  following  relationship 

fi+fv=3fN  (1) 

within  the  accuracy  of  ±1.5 Hz .  Hence,  it  is  quite  possible  that  either  one  of  i  or  i  components  is 
substantial  and  the  other  appears  as  a  result  of  a  nonlinear  interaction  with  the  blade  passing  component. 


Figure  3(b)  shows  the  phase  difference  Arg(pe  / p0)  plotted  against  the  angular  separation  6  of  the 
pressure  taps  (see  Fig.2),  where  the  pressure  fluctuation  is  represented  by  p0  =  pe  exp(2;r//0  .  By 
definition,  the  total  amount  of  the  continuous  change  of  Arg(p0  /  p0)  corresponding  to  the  change  in  6  of 
271,  is  2n  times  the  number  of  cell  n.  If  Arg(pe  /  p0)  decreases/increases  with  an  increase  in  6 ,  it  shows 

that  the  pressure  pattern  rotates  in  the  direction  of/opposite  to  the  impeller  rotation.  We  represent  the 
direction  by  the  plus/minus  signs  on  n  for  co-/counter  rotating  patterns.  From  the  plots  similar  to  Fig. 3 
(b),  the  value  of  n'  for  the  i  component  is  also  determined.  Then  we  can  obtain  the  following  relationship 
consistent  to  the  results  of  nonlinear  interaction  considerations. 

n  +  n'=Zi=3  (2) 

The  propagation  velocity  ratio,  defined  as  the  ratio  of  the  rotational  velocity  of  pressure  pattern  to  the 
impeller  rotational  frequency  is  given  by 

fKnfN)  (3) 
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i  Cavitation  in  backflow  vortices 
(n=5>f/DfN=0-l6) 


it  Rotating  or vita  twn 

(n=Lf/afK=0.9) 


iii  Attached  uorveu  cavitation 

(B=hffnf^U» 


n  Rotating  cavitation 

(tiwltfnf^l.2) 


v  Cavitation  in  backflow  vortices 

(n*5,ffnfN^0.2I) 


Figure  3  (b):  Phase  Difference  Arg(pe  I  p0)  Plotted  against 
the  Angular  Separation  6  of  the  Pressure  Taps. 

The  number  of  cell  n  and  the  propagation  velocity  ratio  f  /(nfN)  thus  determined  are  shown  in  Fig. 3(b). 
We  cannot  distinguish  the  “substantial”  component  from  the  pressure  measurements  alone.  However, 
from  the  flow  visualization  and  other  studies  described  later,  the  components  i  -  v  appear  to  be  the 
substantial  components. 

1.3  Map  of  Oscillating  Cavitation 

Figure  4  presents  the  cascade  plot  of  inlet  pressure  power  spectra  for  various  values  of  static  pressure 
coefficient  y/s .  The  design  point  is  y/s  =  0.130 .  The  number  of  cells  and  the  propagation  velocity  ratio  for 
each  component  are  summarized  in  Table  1.  The  ranges  of  occurrence  of  each  component  are  shown  in 
the  suction  performance  plots  of  Fig. 5.  We  should  note  that  all  of  the  instabilities  occur  in  the  region  of 
cavitation  number  where  the  head  is  not  decreased  by  cavitation.  In  addition,  they  occur  both  above  and 
below  design  flow  coefficient.  These  are  important  characteristics  of  cavitation  instabilities. 


8-4 


RTO-EN-AVT-1 43 


E 

NATO 

3 

OTAN 

Cavitation  Instabilities  in  Inducers 


*  ,  -Q<1£3 


0  I  “0,060 


Figure  4:  Spectra  of  Inlet  Pressure  Fluctuations  for  Various  Static  Pressure  Rise  Coefficient. 


Table  1 :  Number  of  Cells  and  Propagation  Velocity  Ratio 


Number  of  cells 

n 

Rotational  velocity  ratio 
f/nf  N 

i 

Cavitation  in  backflow 
vortices 

5 

0.16  ! 

•  • 

li 

Rotating  cavitation 

1 

*0.9 

... 

in 

Attached  uneven 
cavitation 

1 

1.0 

iv 

Rotating  cavitation 

1 

1. 1-1.3 

V 

Cavitation  in  backflow 
vortices 

5 

0.21  I  j 

vi 

Surge  mode  oscillation 

0 

*18Hz/  (3300-5000  rpm)  j 

vii 

Cavitation  surge 

0 

8— 19Hz/4000  rpm, 
proportional  to  rpm 
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Figure  5:  Suction  Performance  and  a  Map  of  Various  Oscillating  Cavitation  Types. 


1.4  Interpretation  of  Oscillating  Cavitation 
Cavitation  in  Backflow  Vortices. 

Figure  6(a)  presents  a  picture  from  high-speed  video  under  the  condition  with  the  component  i  at 
y/s  =  0.130  and  cr= 0.015.  We  observe  five  clouds  of  cavitation  extending  upstream  from  the  inducer  inlet. 
They  are  supposed  to  be  formed  in  the  backflow  vortices  as  observed  by  Acosta  [2].  The  cavitation  clouds 
slowly  rotate  at  an  angular  velocity  close  to  that  corresponding  to  / /( nfN)  =  0.16 .  Hence,  the  component 
i  is  considered  to  be  caused  by  the  backflow  vortices.  Figure  6(b)  shows  the  picture  at  a  larger  cavitation 
number  <t=0.07  where  the  component  v  is  found.  We  find  a  system  of  cavitation  in  the  tip  leakage  flow, 
which  basically  rotates  attached  to  the  blades,  and  that  in  backflow  vortices  formed  on  the  boundary  of  the 
backflow  region,  rotating  much  slower  than  the  impeller,  as  sketched  in  Fig.  6(c).  As  shown  in  Fig  6(b), 
the  backflow  vortices  occur  rather  irregularly  and  it  was  difficult  to  determine  the  its  number  and  speed 
definitely  from  the  video  picture.  However,  they  were  found  to  be  consistent  with  the  values  n  =  5  and 
//(^/)v)  =  0-2l  determined  from  the  pressure  measurements.  The  passage  of  the  backflow  vortices  is 
considered  to  be  the  cause  of  the  component  v .  Although  the  structure  similar  to  that  as  sketched  in  Fig. 
6(c)  is  observed  at  other  conditions,  distinct  pressure  fluctuation  could  be  found  only  in  the  regions  of  i 
and  v  shown  in  Fig. 5. 
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Figure  6(a):  Cavitation  in  Backflow  Vortices, 
Component  i ,  at  Smaller  Flow  Rate. 


Figure  6(b):  Cavitation  in  Backflow  Vortices, 
Component  v ,  at  Larger  Flow  Rate. 


Tip  leakage  flow  cavitation 


Backflow  region  with  high  tangential  velocity 


Figure  6(c):  Tip  Leakage  Flow  Cavitation  and  Backflow  Vortex  Cavitation. 


Forward  Travelling  Rotating  Cavitation. 

Figure  7  shows  the  fluctuation  of  cavity  length  /,  on  three  blades  of  the  impeller  under  the  condition  with 
the  component  iv .  Forward  propagation  of  cavity  pattern  is  clearly  shown  and  hence  the  component  iv  is 
the  normal  rotating  cavitation.  It  should  be  noted  that  the  maximum  cavity  length  is  slightly  larger  than 
the  circumferential  blade  spacing  H  and  the  mean  value  of  lc!H  is  about  0.75.  As  shown  in  Fig. 5, 
rotating  cavitation  iv  appears  in  two  separate  regions  with  y/s  larger/smaller  than  the  design  value.  In  the 
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region  with  higher  static  pressure  coefficient,  the  length  of  blade  cavity  is  smaller  and  the  cavitating 
region  extends  more  upstream.  With  the  attached  cavitation  component  Hi ,  we  observe  one  shorter  cavity 
and  two  longer  cavities  fixed  to  the  three  blades  of  the  impeller.  Hence,  the  component  Hi  is  caused  by 
unequal  cavities  attached  to  the  rotor.  With  the  component  ii ,  the  change  in  the  size  of  the  cavities  is  not 
so  clear  as  with  the  components  Hi  and  iv .  However,  the  phase  plots  in  Fig. 3(b)  for  a  is  quite  similar  to 
that  for  Hi  and  iv .  This  resemblance  suggests  that  the  component  ii  with  n  =  1  and  f/(nfN)  =  0.9  is 

substantial.  If  so,  this  could  be  the  “backward  traveling  mode”  predicted  by  the  1-D  stability  analysis  [3]- 
[5]  presented  in  another  lecture  of  this  series. 


Figure  7:  Oscillation  of  Cavity  Length  under  Rotating  Cavitation 
(Component  iv ,  y/s  =  0.08  ,  cr  =  0.041 ). 


Surge  Mode  Oscillation. 

For  the  components  vi ,  vii ,  and  vii'  ,  the  pressure  fluctuations  are  in  phase  at  all  circumferential  locations. 
They  are  called  herein  “surge  mode  oscillation”  after  “surge”  in  gas  handling  turbomachines.  The 
fluctuation  of  the  cavity  length  is  plotted  in  Fig. 8  under  the  surge  mode  oscillation  vi .  The  major 
differences  between  the  component  vi  and  the  components  vii  and  vii  are  that,  for  the  component  vi ,  the 
amplitude  is  larger  and  the  frequency  is  fairly  constant  (18  Hz)  as  shown  in  Fig.4.  Figure  9  shows  the 
effects  of  rotational  frequency  fN  at  y/s  =  0.08 .  A  detailed  examination  shows  that  the  frequencies  of  iv 
and  vii  are  proportional  to  fN ,  while  that  of  vi  is  fixed  at  18  Hz.  The  proportionality  of  frequeny  to  the 
rotational  speed  is  an  important  characteristic  of  cavitation  instabilities.  At  y/s  =  0.08  ,  it  was  also  found 
that  fiv  -  fN  is  equal  to  fviV  within  the  measurement  error  (±1.5//z).  Hence,  the  component  vii  is  caused 
by  a  nonlinear  interaction  of  component  iv  (rotating  cavitation)  with  the  rotational  frequency  component. 
The  difference,  fiv  -  fN ,  corresponds  to  the  frequency  of  cavity  oscillation  observed  on  a  blade.  Figure  9 
shows  that  the  component  vi  occurs  when  the  frequency  of  vii '  approaches  18  Hz.  This  suggests  that 
component  vi  is  a  result  of  the  resonance  of  a  certain  vibration  mode  with  18  Hz  and  the  cavity  oscillation 
due  to  rotating  cavitation.  Various  attempts  have  been  made  to  identify  the  “vibration  mode,”  but  no  clear 
mode  could  be  found. 
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ROTATION 


Figure  8:  Oscillation  of  Cavity  Length,  under  Surge  Mode 
Oscillation  (Component  vi ,  y/s  =  0.08  ,  <j  =  0.054  ). 


^S=0.080 


Figure  9:  Effects  of  Rotational  Speed  N  on  the  Frequencies  of  Rotating 
Cavitation  iv  ,  Surge  Mode  Oscillation  vi ,  and  Cavitation  Surge  vii' . 


Although  the  frequency  of  vii  is  also  close  to  fiv  -  fN ,  a  meaningful  difference  (larger  than  ±1.5 Hz ,  up  to 
±5 Hz )  was  observed.  In  addition,  the  component  vii  can  appear  even  without  the  rotating  cavitation 
component  iv .  Figure  10  shows  the  effects  of  the  rotational  speed,  N,  on  the  frequency  of  the  component 
vii .  The  frequency  is  nearly  proportional  to  N  at  lower  frequencies  but  the  rate  of  increase  decreases  at 


RTO-EN-AVT-1 43 


8-9 


Cavitation  Instabilities  in  Inducers 


higher  frequencies.  The  component  vii  is  considered  to  be  a  substantial  one  and  herein  called  “cavitation 
surge”,  since  the  proportionality  of  the  frequency  to  the  rotational  speed  is  an  important  characteristics  of 
cavitation  surge.  No  “resonance”  (or,  a  component  equivalent  to  vi)  was  observed  even  when  fvii 
approached  18  Hz.  This  suggests  that  the  component  vi  is  not  a  simple  “resonance”. 


^S=0.136 

+  a  =0.045 
▼  a  =0.055 
O  a  =0.065 
V  a  =0.075 
▲  o  =0.085 
♦  a  =0.095 
A  a  =0.1 05 
■  a  =0.1 20 
□  o  =0.140 
x  a  =0.150 


0  2000  4000 

N(rpm) 


Figure  10:  Effects  of  Rotational  Speed  N,  on  the  Frequencies  of  Cavitation  Surge  vii . 
The  dash  dot  lines  show  proportionality  relations  between  N  and  f. 


The  modes  of  pressure  fluctuations  iv ,  vi,  and  vii  in  the  inlet  and  outlet  pipes  are  shown  in  Fig.  11.  For  the 
surge  mode  oscillations  vi  and  vii ,  the  phase  is  nearly  constant  throughout  the  pipes  and  the  amplitude 
decreases  linearly  as  the  tank  is  approached  from  the  inducer.  The  pressure  fluctuation  amplitude  of  the 
rotating  cavitation  component  iv  is  significantly  smaller  at  the  inducer  outlet.  It  is  quite  surprising  that  the 
rotating  cavitation  component  iv  can  be  observed  as  far  downstream  as  33  impeller  diameters  from  the 
impeller. 
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Flow  valve  Inducer 


(a)  Rotating  Cavitation  iv  (y/s  =0.080, 
cr  =  0.035 ,  Point  (a)  in  Fig.4). 


(b)  Surge  Mode  Oscillation  vi 

( y/s  =  0.080 ,  cr  =  0.060  5  Point  (b)  in  Fig.4). 


Flow 

<r= 


(c)  Cavitation  Surge  vii  ( y/s  =  0.165  ,  cr  =  0.087  ,  Point  (c)  in  Fig.4). 

Figure  11:  Modes  of  Pressure  Fluctuations  iv,vi ,  and  vii  Represented  by 
the  Amplitude  and  the  Phase  at  Each  Axial  Locations  Shown  in  Fig.1. 


1.5  Effects  of  Piping  System 

In  order  to  examine  the  effects  of  piping  system,  experiments  were  carried  out  with  the  following  two 
additional  configurations  (see  Fig.l): 

•  Diameter  of  the  inlet  pipe  A  is  decreased  from  200mm  to  150mm  (Case  I). 

•  An  accumulator  B  filled  with  air  is  added  (Case  II). 

The  frequency  of  each  component  is  plotted  in  Fig.  12.  The  surge  mode  oscillation  vi  disappeared  when 
we  added  the  accumulator  B.  We  observe  no  significant  change  in  the  frequencies  except  for  the 
cavitation  surge  vii .  This  shows  that  all  of  the  rotating  cavitation  types  i-v  are  quite  independent  on  the 
system.  The  frequency  of  cavitation  surge  vii  is  decreased  by  the  decrease  of  the  inlet  pipe  diameter 
(increase  of  effective  inlet  pipe  length)  and  increased  by  the  addition  of  the  tank  B  (decrease  of  effective 
outlet  pipe  length).  Although  the  frequency  of  cavitation  surge  is  system  dependent  as  shown  above, 
Brennen  [6]  proposes  an  empirical  relation  between  f  /  fN  and  a  for  cavitation  surge.  This  relation  is 
shown  by  the  dashed  curve  in  Fig.  12.  The  present  result  agrees  qualitatively  with  the  empirical  relation. 
The  frequency  of  the  surge  mode  oscillation  vi  was  not  altered  by  the  change  of  the  inlet  pipe  but  it 
disappeared  if  we  added  the  accumulator  to  the  outlet  line.  It  disappeared  also  when  a  small  amount  of  air 
was  introduced  to  the  outlet  line.  These  observations  suggest  that  the  component  vi  is  closely  related  with 
the  outlet  line  but  it  is  not  fully  understood  yet. 
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Figure  12:  Effect  of  the  Piping  System  on  the  Frequencies  of  Oscillating  Cavitation. 
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2.0  CAVITATION  INSTABILITY  CAUSED  BY  PERFORMANCE 
DEGRADATION  DUE  TO  CAVITATION  -  ROTATING  CHOKE 

Most  of  cavitation  instabilities  occur  in  a  range  where  the  performance  degradation  due  to  cavitation  is  not 
significant,  as  shown  in  the  last  section.  However,  in  a  certain  circumstance  near  head  breakdown, 
another  type  of  cavitation  instability  caused  by  the  positive  slope  of  the  performance  can  occur.  This  is 
treated  in  this  section. 

2.1  Experimental  Observation  of  Rotating  Choke 

During  the  firing  test  of  the  LE-7A  LOX/LH2  engine  of  the  HIIA  rocket,  a  large  amplitude  rotor  vibration 
of  the  fuel  turbopump  occurred  at  a  frequency  of  about  350Hz  which  is  about  one-half  of  the  shaft 
rotational  speed  Q  =700Hz.  The  vibration  occurred  when  the  pump  inlet  pressure  was  reduced.  This 
vibration  caused  the  failure  of  the  bolts  fastening  the  bearing  cartridge.  To  investigate  into  this 
phenomenon,  a  series  of  tests  was  conducted  at  full  speed  with  liquid  hydrogen.  The  test  results  are 
reported  by  Shimura  et  al.[7]. 

Figure  13  shows  the  suction  performance  and  the  amplitude  of  the  shaft  vibration  for  two  groups  of  flow 
rates  around  Q/Qd=  0.98  and  Q/Qd=  0.95  .  For  Q/Qd=0.9 5  ,  the  amplitude  becomes  larger  at  the 
cavitation  number  where  the  inducer  head  decreases  rapidly.  The  head  is  kept  nearly  constant,  above  and 
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below  the  cavitation  number.  The  head  decrease  is  not  caused  by  the  instability  since  the  head  decrease  is 
not  recovered  at  smaller  cavitation  number  where  the  amplitude  of  vibration  is  smaller  again.  At  a 
cavitation  number  higher  than  the  breakdown  cavitation  number,  crhigh ,  the  head  is  larger  for  smaller  flow 

rate  and  the  performance  curve  has  a  negative  slope  in  the  flow  rate  vs.  head  plane.  On  the  other  hand,  at 
a  lower  cavitation  number  than  the  breakdown,  cr/ow,  the  head  is  smaller  for  smaller  flow  rate  and  the 
performance  curve  has  a  positive  slope.  These  relations  are  sketched  in  Fig.  14.  For  compressors  and  fans, 
it  is  well  known  that  the  positive  slope  of  the  performance  curve  can  cause  surge  and  rotating  stall  at  a 
smaller  flow  rate. 


Figure  13:  Suction  Performance  and  Shaft  Vibration  with  Two  Groups  of  Flow  Rate. 
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Without  cavitation,  inducers  with  smaller  blade  angle  from  tangent  never  stall  and  the  performance  curve 
has  positive  slope  for  all  flow  rate.  At  higher  cavitation  number  crhigh  where  the  cavity  does  not  extend 

into  the  flow  passage,  no  head  decrease  occurs.  The  head  starts  to  decrease  when  the  cavity  extends  into 
the  flow  channel  between  blades.  It  was  shown  by  Stripling  and  Acosta  [8]  that  the  head  decrease  due  to 
cavitation  can  be  explained  by  the  mixing  loss  downstream  of  the  cavity  terminus.  This  is  called  “choke 
(by  cavitation)”.  The  cavity  thickness  is  larger  for  smaller  flow  rate  and  hence  the  head  decrease  due  to 
cavitation  is  larger,  as  sketched  in  Fig.  15.  If  the  head  decrease  due  to  cavitation  associated  with  the 
reduction  of  flow  rate  is  larger  than  the  increase  of  Euler’s  head,  the  performance  curve  will  have  a 
positive  slope.  So,  it  is  quite  possible  that  the  positive  slope  at  smaller  cavitation  number  is  caused  by  the 
choke. 


Ohigh  c,  Low  Q/Qd 

A,  High  Q/Qd 


Choke 


Glow 


D,  Low  Q/Qd 


B,  High  Q/Qd 


Figure  15:  Sketch  of  Cavity  Shape  at  Higher  and  Lower  Cavitation  Number. 

Figure  16  shows  the  spectrum  of  the  inlet  pressure  fluctuation  and  the  phase  difference  of  the  signals  from 
two  pressure  transducers  located  apart  by  144  degrees  circumferentially  (Shimura,  et  al.,  [7]).  The 
rotational  frequency  of  the  impeller  is  723Hz.  The  phase  difference  at  the  spectrum  peak,  366Hz,  is  147 
degrees,  which  is  close  to  the  geometrical  separation  of  144  degrees.  This  suggests  that  a  disturbance  with 
one  cell  is  rotating  around  the  rotor  at  366Hz,  with  the  rotational  speed  ratio  366/723=0.506.  Since  the 
head  decrease  due  to  cavitation  is  caused  by  choke,  it  is  appropriate  to  call  the  instability  “rotating  choke”. 
All  reported  cavitation  instabilities  occur  in  a  range  of  cavitation  number  where  degradation  of 
performance  due  to  cavitation  is  limited.  So,  this  is  a  new  type  of  cavitation  instability  caused  by  the 
positive  slope  of  pressure  performance  due  to  cavitation  choke.  We  should  note  here  that  the  positive 
slope  of  suction  performance  has  a  stabilizing  effect,  as  shown  by  Young  et  al.  [9]  for  surge  mode 
oscillation. 


Figure  16:  Spectrum  of  Inlet  Pressure  Fluctuation  and  the 
Phase  Difference  between  Two  Pressure  Sensors. 
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2.2  Theoretical  Analysis  of  Rotating  Choke 

Since  the  instability  is  associated  with  the  head  breakdown  due  to  cavitation,  the  analytical  model  should 
be  able  to  simulate  the  head  breakdown.  The  closed  cavity  model  cannot  predict  the  head  decrease  due  to 
cavitation  and  the  results  of  stability  analysis  shown  in  Figure  9  of  the  companion  lecture  [5]  does  not 
include  such  instabilities  associated  with  choke.  To  simulate  rotating  choke,  a  cavity  model  with  a  cavity 
wake  is  applied  as  shown  in  Fig.  17. 


Figure  17:  Waked  Cavity  Model. 

The  viscous/inviscid  interaction  approach  developed  for  cavitation-free  flows  is  used.  According  to  this 
approach  the  flow  is  divided  into  an  internal  turbulent  flow  behind  the  cavity  and  an  external  inviscid  flow 
in  which  the  cavity  is  closed  on  the  body  of  viscous  layer  displacement  (Semenov  et  al.,  [10]).  In  the 
cavity  wake  three  specific  regions  are  considered.  The  first  is  the  transitional  mixing  region  where  the 
velocity  profile  across  the  wake  is  under  formation  and  the  density  changes  smoothly  due  to  flow  re¬ 
circulation.  This  region  is  needed  to  provide  a  good  agreement  of  suction  performance  with  experiment 
(Semenov  et  al.,  [11]).  The  second  is  the  near  wake  in  which  the  velocity  profile  changes  according  to  the 
external  pressure  gradient  along  the  wake.  Von  Karman’s  integral  equation  is  used  and  it  is  assumed  that 
the  coefficient  of  turbulent  mixing  is  constant  (Gogish  et  al.,  [12]).  The  third  is  the  far  wake  which  starts 
from  the  trailing  edge  of  blades  and  adjusts  itself  on  the  downstream  flow.  The  conditions  of  interaction 
between  viscous  and  inviscid  flows  in  the  near  wake  make  it  possible  to  find  the  shape  of  the  boundary  of 
the  inviscid  flow  and  to  calculate  both  steady  and  unsteady  flows.  The  normal  velocity  of  the  outer  flow 
on  the  boundary  is  equated  with  that  determined  from  the  change  of  the  wake  thickness. 

The  head  decrease  due  to  cavitation  is  associated  with  the  increase  of  relative  velocity  caused  by  the 
displacement  effect  of  the  cavity  wake.  The  stability  analysis  is  similar  to  that  with  closed  cavity  model 
(Horiguchi  et  al.,  [13], [14])  described  in  [5].  The  problems  for  the  inviscid  and  viscous  flow  are  reduced  to 
a  system  of  linear  equations  by  using  linear  interpolation  of  source  and  vortex  and  finite  difference 
approximation  of  the  ordinary  differential  equations  for  the  viscous  wake.  The  solution  of  the  eigenvalue 
problem  makes  it  possible  to  determine  the  frequencies  and  the  stability  of  the  various  modes  of  cavitation 
instabilities. 

Figure  18  shows  the  suction  performance  predicted  by  the  model.  On  the  suction  performance  curves,  the 
occurrence  of  rotating  cavitation  is  shown  by  open  symbols,  while  the  occurrence  of  rotating  choke  is 
shown  by  closed  symbols.  Other  types  of  instabilities  such  as  shown  in  Fig. 9  of  [5]  are  also  predicted  by 
the  present  model  but  they  are  not  shown  here.  We  note  that  the  rotating  cavitation  mainly  occurs  in  the 
range  where  the  performance  degradation  is  insignificant,  while  rotating  choke  occurs  where  the  head  is 
decreased  due  to  cavitation. 
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Figure  18:  Suction  Performance  Predicted  by  the  Waked  Cavity  Model. 


The  performance  curve  in  the  flow  coefficient  vs.  head  plane  is  shown  in  Fig.  19,  where  the  occurrences  of 
rotating  cavitation  and  rotating  choke  are  also  shown  by  open  and  closed  symbols,  respectively.  The  same 
observations  as  with  Fig.  18  are  also  evident  in  this  figure.  Rotating  choke  occurs  typically  in  the  region 
with  positive  slope  of  the  performance  curve,  which  is  quite  different  from  rotating  cavitation.  This  clearly 
shows  that  rotating  choke  is  caused  by  the  positive  slope  of  the  performance  curve  due  to  choke. 
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Figure  19:  Head  Coefficient-Flow  Coefficient  Curves  for  Various  Cavitation  Number. 


Figure  20  shows  the  propagation  velocity  ratio  of  rotating  cavitation  and  rotating  choke.  Rotating  choke  is 
found  only  at  a  region  with  smaller  cj  .  Rotating  cavitation  rotates  faster  than  impeller  but  rotating  choke 
rotates  about  a  half  of  the  impeller  speed.  Rotating  stall  in  impellers  with  longer  blade  rotates  with  a 
speed  only  slightly  smaller  than  the  impeller  speed,  caused  by  the  inertia  effect  of  the  fluids  in  the  impeller 
(Tsujimoto,  et  al.,  [3], [4]).  Lower  speed  of  the  rotating  choke  may  be  caused  by  the  decrease  of  the  inertia 
effect  due  to  cavitation. 
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Figure  20:  Propagation  Velocity  Ratio  of  Rotating  Cavitation  and  Rotating  Choke. 

Figure  21  shows  the  cavity  shapes  under  rotating  cavitation  and  rotating  choke  at  their  onset  points  (the 
largest  cavitation  number).  For  rotating  choke  the  mixing  region  extends  to  the  throat  of  the  blade  channel. 
On  the  other  hand,  rotating  cavitation  starts  to  occur  with  much  shorter  cavity.  With  a  closed  cavity  model, 
it  has  been  shown  that  various  types  of  cavitation  instability  start  to  occur  when  the  cavity  length  reaches 
about  65%  of  the  spacing  as  mentioned  in  [5].  Figure  22  shows  the  steady  cavity  length  lciosecj  of  the 
closed  cavity  model  and  the  total  length  of  cavity  and  mixing  region  lc  +  lm  of  the  waked  cavity  model, 
normalized  by  the  blade  spacing  h ,  plotted  against  crlila) .  The  regions  of  rotating  cavitation  and  rotating 
choke  occurrence  are  also  shown.  It  is  clearly  shown  that  rotating  cavitation  starts  to  occur  when  the 
cavity  length  (closed  cavity  model)  or  the  total  length  of  cavity  and  mixing  region  (waked  cavity  model) 
exceeds  65%  of  the  blade  spacing,  h .  This  suggests  that  the  total  length  of  cavity  plus  mixing  region  plays 
an  important  role  in  the  present  model.  On  the  other  hand,  rotating  choke  starts  to  occur  when  the  total 
length  of  cavity  and  mixing  region  exceeds  150%  of  the  spacing. 


a)  rotating  cavitation  (cr  =  0.26)  b)  rotating  choke  (a  =  0.077) 

Figure  21 :  Cavity  Shapes  at  the  Onset  Point  of  Rotating  Cavitation  and  Rotating  Choke  at 
(j)  =  0.095 .  Cavity  is  shown  by  the  white  area.  The  density  in  the  mixing  region  and  the 
wake  is  shown  by  the  darkness.  Vertical  dimension  is  magnified  with  a  factor  3. 
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Figure  22:  Cavity  Length  of  Closed  Cavity  Model  lciosecj ,  and  the  Total  Length  of  Cavity  and 
Mixing  Region  lc  +  lm  of  the  Waked  Cavity  Model.  The  gray  bars  show  the  range  of 

rotating  cavitation  and  rotating  choke.  (J3  =  80°  ,  C/h  =  2.35  ,a  =  4°). 


Thus,  the  present  model  can  simulate  two  important  characteristics  of  rotating  choke,  the  onset  condition 
and  frequency.  However,  predicted  suction  performance  shown  in  Fig.  18  is  significantly  different  from  the 
experimental  curve  as  shown  in  Fig.  13.  This  may  be  caused  by  non-linear  and/or  three-dimensional 
effects.  Positive  slope  and  hence  rotating  choke  are  unusual  and  not  observed  for  many  inducers.  Further 
research  is  needed  for  the  complete  understanding  of  rotating  choke. 


3.0  HIGHER  ORDER  MODES  OF  CAVITATION  INSTABILITIES 

In  addition  to  cavitation  surge  and  forward  mode  of  rotating  cavitation,  various  kinds  of  cavitation 
instabilities  are  predicted  by  a  2-dimensional  flow  stability  analysis  of  a  cavitating  cascade  with  a  closed 
cavity  model  [5].  These  modes  were  subsequently  found  by  experiments  and  they  are  treated  in  this 
section  (Tsujimoto  et  al.,  [15]). 

3.1  Backward  Rotating  Cavitation 

Backward  travelling  rotating  cavitation  was  identified  for  the  first  time  by  Hashimoto  et  al.  [16].  Figure 
23  shows  the  spectrum  of  inlet  pressure  fluctuation  of  a  3-bladed  inducer  rotating  at  fN  =  \\lHz. 
Component  F  with  fF  =  13SHz  is  conventional  forward  propagating  rotating  cavitation.  With  this 
component,  an  additional  component  with  3 (fF  -  fN)=64Hz ,  which  is  the  frequency  at  which  the  blades  cut 
the  cavitating  region,  is  observed.  At  the  cavitation  number  slightly  larger  than  the  forward  rotating 
cavitation  onset,  Component  B  with  fB  =  159 Hz  is  observed.  With  this  component,  an  additional 
component  with  3 (fB  +  fN)  =  828/fz  is  observed.  This  frequency  is  obtained  if  we  assume  that  the 
cavitated  region  rotates  in  the  direction  opposite  to  the  impeller  rotation.  Observation  of  cavity  oscillation 
by  high-speed  film  also  suggests  backward  rotation. 
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Figure  23:  Spectrum  of  Inlet  Pressure  Fluctuation  of  a  3-Bladed 
Inducer  Showing  Forward  and  Backward  Rotating  Cavitation. 

Another  example  of  backward  rotating  cavitation  was  obtained  during  the  tests  of  a  forward  swept  4- 
bladed  inducer.  Figure  24  (a)  shows  the  spectrum  of  inlet  pressure  fluctuation.  The  rotational  frequency  is 
fN  =  66.1Hz  .  We  observe  conventional  rotating  cavitation  around  94Hz  (/  lfN  =  1.41)  for  the  entire  range 
of  cavitation  number  shown.  Alternate  blade  cavitation  occurs  around  cr=  0.063.  Backward  rotating 
cavitation  is  found  at  129Hz  (/  lfN  =  (— )1 .95 )  for  0.065  <cr<  0.078.  This  propagation  velocity  ratio  of  - 
1.95  is  larger  than  that  for  the  three  bladed  inducer  (-159/1 17=-1. 36)  mentioned  before,  and  is  in 
qualitative  agreement  with  the  theoretical  results  of  -1.25  for  four  bladed  inducer  and  -1.0  for  three 
bladed  inducer  at  a /2a  =  1.0 . 
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Figure  24:  Backward  Rotating  Cavitation  Observed  in  a  Forward  Swept  4-Bladed  Inducer. 

Figure  24  (b)  shows  the  plot  of  the  phase  of  inlet  pressure  fluctuation  components  at  various 
circumferential  locations  denoted  by  6.  For  the  backward  rotating  cavitation,  the  phase  advances  as  we 
proceed  in  the  direction  of  the  impeller  rotation  and  reaches  360  degrees  as  we  make  one  turn.  This  shows 
that  a  pressure  pattern  with  one  cell  is  rotating  in  the  direction  opposite  to  the  impeller.  Fig.  24(b)  also 
shows  that  the  conventional  rotating  cavitation  and  the  alternate  blade  cavitation  have  one  and  two  cells 
respectively  and  that  both  of  them  are  rotating  in  the  same  direction  as  the  impeller. 

3.2  Higher  Order  Rotating  Cavitation 

Figure  25  (a)  shows  the  spectrum  of  the  inlet  pressure  fluctuation  observed  for  a  3-bladed  inducer  at  the 
rotational  speed  fN  =  50Hz  .  Conventional  rotating  cavitation  is  observed  for  0.04  <  cr<  0.058  ,  at 
fc  =  61.25 Hz .  In  the  region  of  cavitation  number  larger  than  that  of  conventional  rotating  cavitation  onset, 
0.06  <  <j  <  0.09  ,  we  observe  a  component  with  fh  o  =  250 Hz  .  Although  hidden  behind  the  250Hz 
component,  there  is  a  component  with  fhs  =  243/ 15Hz  over  0.058  <  <j  <  0.068  .  Figure  25(b)  shows  the 
spectrum  of  stress  fluctuation  of  a  blade.  The  components  with  nearly  constant  frequency  are  considered 
to  be  electrical  noise.  With  conventional  rotating  cavitation,  for  0.04  <  cr<  0.058,  we  observe  a  component 
with  fc  —  fN  =  11.25ft  which  is  the  frequency  of  conventional  rotating  cavitation  observed  in  the  rotating 
frame.  In  the  same  way,  we  observe  a  component  with  fho-fN  =  200Hz,  which  is  the  frequency  of 
forward  propagating  fho  =  250 Hz  component  observed  in  the  rotating  frame.  So,  this  component  may 
correspond  to  the  higher  order  rotating  cavitation.  In  the  region  0.06  <  a  <  0.068 ,  we  observe  a  component 
with  fhs  =  243.15Hz  ,  which  is  the  same  as  that  observed  in  the  inlet  pressure  fluctuation.  This  can  be  the 
higher  order  surge  mode  oscillation. 
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(a)  Spectrum  of  inlet  pressure  fluctuation 


(b)  Spectrum  of  stress  fluctuation  of  a  blade 


Figure  25:  Spectrum  of  Inlet  Pressure  and  Blade  Stress  Fluctuations,  Showing  a  Higher  Order 
Rotating  Cavitation  and  a  Higher  Order  Surge  Mode  Oscillation. 


Figures  26  (a)  and  (b)  show  the  phase  (positive  value  means  phase  advance  relative  to  PI  and  Gl)  of  the 
pressure  and  stress  fluctuations  measured  at  various  circumferential  locations.  They  clearly  show  that  both 
conventional  rotating  cavitation  (with  61.25Hz  and  11.25Hz)  and  the  higher  order  rotating  cavitation  (with 
250Hz  and  200Hz)  have  one  cell  and  rotate  in  the  direction  of  impeller  rotation.  The  higher  order  surge 
mode  oscillation  with  fhs  =  243.15Hz  has  the  same  phase  over  the  circumference.  The  frequencies  of 
higher  order  rotating  cavitation  and  higher  order  surge  for  3-bladed  inducer  predicted  by  the  2-D  stability 
analysis  [5]  are  both  about  4.2  times  the  rotational  frequency  at  cr/2a:  =  1.0  .  This  is  close  to  the 
experimental  values  of  5.0  and  4.88  for  higher  order  rotating  cavitation  and  higher  order  surge  mode 
oscillation  respectively. 
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2 

(b)  Blade  stress  fluctuations 


Figure  26:  Phase  of  Inlet  Pressure  and  Blade  Stress  Fluctuations, 
Corresponding  to  the  Spectrum  in  Fig.25. 
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The  component  in  Fig. 25  (a)  around  240Hz  observed  over  cr>  0.075  has  low  coherence  between  the 
pressure  signals  at  different  circumferential  locations.  For  this  reason,  and  also  because  it  also  occurs  at 
higher  inlet  pressure  without  cavitation,  we  consider  that  the  component  is  caused  by  a  backflow  vortex 
structure  at  the  inlet. 

3.3  Higher  Order  Surge  Mode  Oscillation 

During  the  engine  firing  and  water  tests  for  the  assessment  of  the  fatigue  failure  of  a  blade  of  the  liquid 
hydrogen  turbopump  inducer  experienced  in  the  8th  launch  of  HII  rocket,  a  high  frequency  pressure 
fluctuation  was  found.  The  frequency  was  about  4.7  times  the  rotational  frequency  which  is  close  to  the 
first  bending  mode  frequency  of  the  blades.  For  this  reason  extensive  tests  were  carried  out  (Tsujimoto  et 
al.,  [17]). 

Figures  27  (a)  and  (b)  show  the  spectrum  of  pressure  fluctuations  measured  at  two  circumferential 
locations  separated  by  90  degrees,  on  the  casing  wall  and  at  mid-chord  obtained  in  a  water  test.  The 
inducer  has  three  blades  and  the  rotational  frequency  is  fN  =  l50Hz.  We  observe  sharp  spectrum  peaks  at 
the  blade  passing  frequency  3 fN  and  its  harmonics.  In  addition,  we  observe  a  broadband  component 
centered  at  705Hz,  which  corresponds  to  4.7  times  the  rotational  frequency.  Note  that  the  peak  is  higher 
than  the  blade  passing  component.  Figure  27  (c)  and  (d)  show  the  coherence  and  the  phase  of  the  cross 
spectrum  between  those  pressure  signals.  We  observe  high  coherence  around  705Hz,  and  the  phase 
difference  there  is  zero.  This  means  that  the  705Hz  component  is  a  surge  mode  oscillation.  Blade  stress 
fluctuation  has  similar  spectrum  and  has  high  coherence  with  the  case  pressure  fluctuation  around  705Hz. 
This  component  is  important  since  the  frequency  can  be  close  to  that  of  the  first  bending  mode  of  blade 
vibration. 
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(a)  Case  pressure  fluctuation  at  midchord 


0  2kHz 

(b)  Case  pressure  fluctuation  at  midchord, 
90  deg  drcumferen Italy  separated  from  (a) 


shown  in  (a)  and  (b) 


(d)  Phase  of  cross  spestrutn  between 
the  signals  shown  in  (a)  and  (b) 


Figure  27:  Spectrum  of  Pressure  Fluctuation  at  Two  Different  Circumferencial 
Locations,  with  Coherence  and  Phase  Difference  between  the  Signals. 


The  following  characteristics  are  found  in  the  series  of  water  tests  for  the  705Hz  component: 

(1)  The  frequency  does  not  depend  on  the  length  of  inlet  pipe.  So,  it  is  system  independent. 

(2)  The  onset  cavitation  number  ( a  =  0.1 )  is  significantly  larger  than  that  of  rotating  cavitation 
( <j=  0.06). 

(3)  The  frequency  is  proportional  to  the  rotational  frequency. 

The  characteristics  of  (1)  to  (3)  agree  with  those  of  the  higher  order  surge  mode  oscillation  (Mode  IX) 
predicted  by  the  2-D  stability  analysis  in  [5]. 

The  Strougal  number  of  Mode  IX  based  on  cavity  length  is  kept  nearly  constant  (Sf  =  1.3)  and  the 
frequency  changes  from  3  times  the  rotational  speed  at  cr/2a:  =  0.12  to  12  times  at  cr/2a:  =  5.61.  This 
includes  the  experimental  value,  4.7  times  the  rotational  frequency,  but  such  large  dependence  on  alia 
was  not  identified  in  experiments. 
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It  is  unfortunate  that  further  assessment  as  to  closer  correlation  of  the  higher  order  surge  mode  oscillation 
with  the  fatigue  failure  of  the  blade.  However,  it  is  important  to  remember  that  we  do  have  high 
frequency  cavitation  instabilities  which  might  resonate  with  the  structural  vibration  modes  of  blades. 

The  backward  rotating  cavitation,  the  higher  order  rotating  cavitation,  and  the  higher  order  surge  mode 
oscillation  introduced  here  are  all  quite  repeatable  for  each  case.  However,  they  can  be  observed  only 
under  limited  circumstances,  as  compared  with  conventional  cavitation  surge  and  rotating  cavitation.  The 
reason  for  that  is  not  known  as  yet. 


4.0  CONCLUSION 

Experimental  observations  of  various  kinds  of  cavitation  instabilities  are  presented.  Typical  cavitation 
instabilities  occur  in  the  range  of  cavitation  number  where  the  inducer  head  is  not  affected  by  cavitation. 
This  type  of  cavitation  instability  is  caused  by  the  cavity  volume  fluctuation  and  includes  cavitation  surge, 
and  forward  rotating  cavitation.  They  occur  even  at  design  flow  coefficient.  The  frequency  of  cavitation 
surge  is  proportional  to  the  rotor  speed  and  the  cavitated  region  rotates  faster  than  impeller  for  forward 
rotating  cavitation. 

Another  type  of  cavitation  instability  is  caused  by  the  positive  slope  of  the  pressure  performance  curve  due 
to  the  blockage  effect  of  cavitation.  This  type  of  cavitation  instability  occurs  at  smaller  cavitation  number 
than  normal  rotating  cavitation  and  the  cavitated  region  rotates  slower  than  impeller. 

Stability  analyses  of  1-  and  2-dimensional  cavitating  flow  predict  various  modes  of  cavitation  instabilities. 
Usually  only  cavitation  surge  and  forward  rotating  cavitation  are  observed  but  other  modes  of  cavitation 
instabilities  such  as  backward  rotating  cavitation,  higher  order  modes  of  cavitation  surge  and  rotating 
cavitation  are  also  observed  only  under  limited  conditions.  Observed  in  a  frame  fixed  to  the  rotor,  these 
higher  order  modes  have  very  high  frequencies  and  they  may  excite  resonance  of  blades. 

Since  cavitation  instabilities  occur  even  at  design  flow  rate  without  head  breakdown,  they  may  lead  to 
disastrous  results  under  continuous  operation.  So,  it  is  important  to  check  that  they  are  adequately 
suppressed  under  all  operating  conditions. 
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